Description of single-A hypernuclei with relativistic point coupling model 
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We extend the relativistic point coupling model to single-A hypernuclei. For this purpose, we add 
N-A effective contact couplings to the model Lagrangian, and determine the parameters by fitting 
to the experimental data for A binding energies. Our model well reproduces the data over a wide 
range of mass region although some of our interactions yield the reverse ordering of the spin-orbit 
partners from that of nucleons for heavy hypernuclei. The consistency of the interaction with the 
quark model predictions is also discussed. 
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I. INTRODUCTION 

Relativistic mean field (RMF) theory has been suc- 
cessfully applied to both finite nuclei and nuclear matter 
in order to describe their bulk properties [H-13|- Start- 
ing from an effective Lagrangian in which nucleon and 
meson fields are coupled in a covariant manner, single- 
particle Dirac equations for nucleons are derived within 
the mean field approximation. In this model, nucleons 
are treated as Dirac particles moving independently in 
the mean field generated by the mesons. The spin-orbit 
interaction with the correct sign and magnitude naturally 
arises from the relativistic treatment of nucleons. The 
success of the model affirms the meson exchange picture 
of a nucleon-nucleon interaction in nuclei. 

Recently, another class of relativistic model, that is, 
the relativistic point coupling (RPC) model y| 0| has 
also been widely employed 0" 13 1 . This model consists of 
Skyrme-type zero-range interactions and has been found 
to be as capable as meson exchange RMF models of 
reproducing the properties of finite nuclei and nuclear 
matter [f| @, fiol — Tl2l ] . This model has several advantages 
compared to the meson exchange models. First, there is 
no need to solve the Klein-Gordon equations for mesons 
since the mesonic degrees of freedom are all implicit in 
the RPC model. Secondly, the Fock terms can easily be 
introduced by using the Fierz transformation because of 
its zero-range nature [H, Q . Lastly, it is much easier to 
apply the model to beyond-mean-field methods such as 
the generator coordinate method (GCM), and angular 
momentum and particle number projections [TTI . Il2| 

A zero-range type interaction is suitable also for the 
three-dimensional (3D) mesh method for mean-field cal- 
culations in the coordinate space representation 0, [HI ■ 
With this method, an arbitrary deformation of nuclei can 
be efficiently described, and the method has been widely 
used in non-relativistic Skyrme Hartree-Fock (SHF) cal- 
culations together with the imaginary time technique 
[iH Il5j . An extension of this method to the relativistic 
approach is not trivial, however. That is, a naive imag- 
inary time evolution breaks down in relativistic systems 
due to the presence of the Dirac sea 0, [0} ■ This is not 
a numerical, but rather a fundamental problem related 
to the variational principle. Such phenomena have been 



well-known under the name of "variational collapse" in 
the field of relativistic quantum chemistry [l8l - |22j . For 
this reason, a 3D mesh calculation has not yet been car- 
ried out neither with RMF nor RPC. Recently, a few 
prescriptions to avoid the variational collapse have been 
tested in the nuclear physics context [H, [jj} • The pre- 
scriptions were found to work well, at least for simple 
spherical systems, and a 3D mesh calculation may now 
be almost ready to perform. 

In this paper, we extend the relativistic point coupling 
model to hypernuclei. Effects of adding a A particle on 
the shape of nuclei, i.e. the glue-like role, are attracting 
much attention theoretically and experimentally [23M23j . 
The RMF theory with meson exchange model has been 
extended to hypernuclei in order to describe such effects 
by adding ordinary scalar and vector coup lin gs o f a A par- 
ticle to a and u> mesons, respectively [23l |24| . |3fJ-|40( ■ The 
authors of Refs. [3(J, HH fitted the coupling constants in 
the strange sector to the experimental data of A binding 
energies. In these calculations, the tensor couplin g to the 
u) meson, which is predicted by the quark model |32l l33j 
to be much stronger for A than for nucleon, was also in- 
cluded in order to reproduce rather small A spin-orbit 
splittings. Our aim in this paper is to propose a zero- 
range version of those phenomenological RMF model for 
hypernuclei. 

The paper is organized as follows. In Sec. [ill we intro- 
duce the model Lagrangian for a point coupling model 
extended to single-A hypernuclei. In Sec. IIII1 the opti- 
mal parameter set is obtained by fitting to experimental 
single-particle energies for the A particle. We also test 
its predictive power and a consistency with the quark 
model. In Sec. IIV1 we summarize the paper. 



II. THE MODEL 



by 



Our model Lagrangian for single-A hypernuclei is given 
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in which the free and electromagnetic parts are given by 



£frcc = ipN(i$ ~ rn N )i} N + ipA{i$ ~ rn A )^A, (2) 

7 V", (3) 



respectively. Here, V'jvj V'a? an d are the nucleon, 
lambda, and electromagnetic fields, respectively. P" = 
d^A" — d v A^ is the electromagnetic field strength tensor. 
The masses of nucleon and lambda arc denoted by 
and m\, respectively. T3 is the isospin matrix. 

The RPC model for normal nuclei consists of four- 
fermion point couplings £4f, derivative terms £dor, and 
higher order terms £hot 0- Here, is the leading-order 
of zero-range approximation to the meson exchange in- 
teraction. £der simulates the finite ranges of the meson 
exchanges. £hot corresponds to the self couplings of the 
scalar and vector mesons, which introduces a density de- 
pendence into N-N contact couplings. These terms are 
given by 



with 
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and 



1 - - 2 

- 4 TV [{^n1^n){'^ni' j "^n)\ ■ 



(7) 



Notice that the four different spin-isospin vertex struc- 
tures labeled by the subscripts S, V, TS, and TV in the 
coupling constants correspond to a, uj, (5, and p meson 
exchanges, respectively. Thus we can find one-to-one cor- 
respondence of each term to the meson exchange model. 



Noticing that A only couples to the scalar and vector 
mesons, we construct AT- A interaction as 
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where 
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and 



cZ = a { T NA) ttJA<j^ip A )(d^ Nl M. 
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For simplicity, we do not consider the higher order 
term for the ATA coupling, C^ t , in this paper. 
in Eq. (fTT|) simulates the A-oj tensor coupling = 
i^j( , ipA& tJ ' u '4>A)(dvU^)- As we mentioned in the Intro- 
duction, the quark model suggests that the tensor cou- 
pling of A to w meson is much stronger than that of 
nucleon. That is, the quark model yields the ratio of A-w 
tensor-to-vector coupling constants, /aw/^Alj, to be — 1, 
while it yields the corresponding ratio for nucleon to be 
fNui/gNu — —0.09 [HI]. Thus this type of coupling plays 
an important role in hypernuclei. Since this term is pro- 
portional to the derivative of the mean field, it mainly af- 
fects the spin-orbit splittings of A single-particle energies 
}30j . It is expected that the small spin-orbit splittings of 
lambda can be reproduced by tuning the tensor coupling 
dj?^ . We will discuss this point in the next Section. 

Our model presented in this paper is similar to the 
one adopted in Ref. [13] by Finelli et al., which is based 
on the chiral SU(3) dynamics. The A^-A part of their 
model consists of the density dependent contact four 
fermion couplings of the scalar and the vector types, 
as well as the derivative term of the scalar type. A 
part of the four fermion terms effectively describe un- 
resolved short distance physics, while the density depen- 
dence is attributed to in-medium Nambu-Goldstone bo- 
son (two-pion and kaon) exchanges. The coefficients for 
the density-dependence are fixed by the chiral SU(3) per- 
turbation theory. The tensor interaction caused by 2n ex- 
change is treated in the model of Ref. [13] as first-order 
perturbation on the Hartree single-particle energies. In 
contrast to the model of Finelli et al. , we determine all the 
parameters phenomenologically, and thus we consider in 
our model also the derivative coupling of the vector type, 
which is absent in theirs. 

The total energy corresponding to the Lagrangian in 
Eq. (TTJ for a single-A hypernucleus with mass number A 
(i.e., a single A particle with A — 1 nucleons) in the mean 
field (Hartree) and the no-sea approximations is given by 
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where a and (3 are the usual Dirac matrices. Here, we 
have assumed the time reversal invariance of the nuclear 
ground state. The densities appearing in Eq. (TI"2"|) are 
defined as 
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Here ipi is the wave function for the i-th nucleon, and ip\ 
is the wave function for the A particle. 

The relativistic Hartree equations for the nucleons and 
lambda particle are obtained by taking the variation of 
the energy with respect to the wave functions as, 

l^( E -p'J d3r '^)= - <17) 

where ti is a Lagrange multiplier which ensures the nor- 
malization of the single particle wave functions. Varia- 
tion with respect to the nucleon wave function leads to 
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the Hartree equation for nucleons, 

[a-p + Vv + VtvTz + Vc+ {m N +V S + V T st 3 ) fi] ipi 

= tiipi, 
(18) 

with 

Vs = a s ps + PsPs + IsPs + S s Ap s 

+4 m) 4 A) + 4 m) A4 A \ 

Vv = aypv + IvPv + SyApy 

+ S^ApP + a (^) p CA) , (19) 

VtS = OiTSPTS + Sts^Pts, 
Vtv = OiTVPTV + Stv Apxv > 
V c = eA 01 -^, (AA° = -e$ } ), 

while variation with respect to the lambda wave function 
leads to the Hartree equation for the lambda particle: 

[3 -p + U v + U T + (m A + U s ) 0\ ipA = caV'a, (20) 
with 

Us = af^ps + S^Aps, 

Uv = a^py + S^Apy, (21) 

U T = -ia ( T A) (3a- {Vp v ). 

After having solved these Hartree equations self- 
consistently, we obtain the total binding energy as 
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where the center of mass energy Eqm is calculated by 
taking the expectation value of the kinetic energy for the 
center of mass motion with respect to the many-body 
ground state wave function as 



Eqm 
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2[(A- l)m N + m A ] 
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See Appendix for the explicit expression for this term. 
The relation of the point coupling model to the me- 



son exchange model can be made as follows (See also 
Eqs. (6)-(10) in Ref.Q). By eliminating the meson fields 
and expanding the meson propagators to the leading or- 
der, the following approximate relations between the two 
models can be obtained Q: 
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where g's and m's are the baryon-meson coupling con- 
stants and the meson masses, respectively. f Au is the 
A-u> tensor coupling constant. Notice that it has been 
demonstrated that as and ay obtained phenomenolog- 
ically approximately follow these relations @ . If we as- 
sume the naive quark counting ratios g Aa = fffAfo- and 
.9Aoj = fffJVw, together with the quark model prediction 
for the tensor coupling, /a^/saw = — 1> we obtain 
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We will show in the next Section that these expected re- 
lations indeed hold if we include the N-A tensor coupling 
given by Eq. in the Lagrangian. 



III. RESULTS AND DISCUSSION 

With the model described in the previous section, we 
calculate A binding energies defined by the mass differ- 
ence 

m( A ~ 1 Z)+m A -m(£z)=E B ( A - 1 Z)-E B (£z). (27) 

To this end, we assume spherical symmetry, and neglect 
the pairing correlations for simplicity. For the valence 
orbit, we use the filling approximation to determine the 
occupation probability. We set the masses of baryons to 
m N = 938 MeV and m A = 1115.6 MeV. We use the 
parameter set PC-F1 [tJ for the N-N part of interaction 
and fit the five parameters in the N-A part (see Eqs. ([5]), 
(fl~0|) and (fTTj) ) to the experimental data. The data to be 
fitted to are A binding energies for s and p orbitals in 
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Ca, s and d in 5 \V, s, p, d, and / in 
8 ^Y, s, p, d, /, and g in 13 ^La, and s } p, ri, /, and g in 



16 n 



s, p and d in 
208 Pb. These are taken from Refs.l41| 



A 

the spin-orbit splitting for the p orbital of A in A 



In addition, 
\0 [H| is 

included in the fitting procedure. The value deduced in 
Ref.[ii| is 300 keV < e Api/2 -e Ap3/2 < 600 keV, where the 
variation comes from a choice of the interactions. Notice 
that this value is model dependent, and we merely regard 
it as a criterion. The coupling constants in the strange 
sector are determined by performing a least-squares fit 
to the data, that is, by minimizing the quantity 
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Here, Adof is the number of degree of freedom, and O* hcor 
and 0® xpt are theoretical and experimental values of the 
observables, respectively, with the experimental uncer- 
tainties of AO^ xpt . To find the minimum of Xdoi m 
the five dimensional parameter space, we employ an au- 
tomatic search algorithm Oak-ridge and Oxford method 
431. 



The parameter set PCY-S1 so obtained is summarized 
in Table HI Together with the coupling constants, the 
ratios R of the resultant A^-A coupling constants to the 
expected values given in Eq. (|2"6")l , 

R = (resulted value) /(expected value), (29) 

are also shown. These ratios are R = 0.79, 0.96, and 
1.37 for ag~\ ay^ A \ and a^ A \ respectively, and the 
expected values are approximately realized. 

The calculated binding energies of A with this inter- 
action are shown in the upper panel of Fig. [1] One ob- 
serves that the calculated A binding energies agree with 
the experimental values fairly well, although the binding 
energies for 2 A Si and 3 A S are somewhat overestimated. 
The less satisfactory result for these latter nuclei, which 
has been observed also in the previous RMF calculations 
for hypernuclei [3(1 [H, is within expectation, 
as we do not take into account a strong deformation of 
the core nucleus nor the pairing correlation. We have 
confirmed that the situation does not change even if we 
include these two nuclei in the fitting. 

In order to investigate the role of the tensor coupling, 
we show in Table [TTI the parameter set PCY-S2 obtained 
without including the tensor coupling term. The Lambda 
binding energies calculated with this interaction is shown 
in the upper panel of Fig. [5] As one sees, the agreement 
with the experimental data is worsened as compared to 
PCY-S1, and the ratios R are strongly suppressed com- 
pared to unity. On the other hand, the sum ct s NA ^ +a{f A ' 
has similar values around -3 x 10" 5 MeV~ 2 for PCY-S1 
and PCY-S2. The suppression of the ratios can be un- 
derstood as follows. In the non-relativistic reduction of a 
Dirac equation without the tensor coupling contribution, 
the central potential and the spin-orbit potential read 

Ventral = V + S, V ls = - ( V - S), (30) 

2m z r dr 

where V and S are the vector and the scalar potentials, 
respectively. Therefore, to reproduce a small spin-orbit 
splitting of A without the tensor interaction, the differ- 
ence of the vector and the scalar potential have to be 
small, keeping their sum constant. This can be achieved 
only by lowering the values of the four fermion A-A cou- 
plings, a s NA ^ and ay^ A \ which roughly determine the 
strengths of mean potential felt by A. Notice that V — S 
does not have to be small in the presence of the tensor 
coupling, as there is another contribution to the spin- 
orbit potential from the tensor coupling. The importance 
of the A^-A tensor coupling (originated from the A-tu ten- 
sor coupling) is thus evident. It yields small spin-orbit 
splittings, keeping a^ A ^ and a^ 1 ^ at the natural values. 
In PCY-S1, the two quark model predictions, that is, the 
quark counting ratios and the importance of the tensor 
coupling (/aw/sAw = — 1), are simultaneously satisfied. 

Let us now discuss the calculated spin-orbit splittings, 
AE so . These are estimated as a difference of A single- 
particle energies between spin-orbit partners, A_B so = 
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TABLE I. The best fit parameter set PCY-S1 for the relativis- 
tic point coupling model for hypernuclei. PC-F1 [7J is used 
for the N-N part. The ratios R defined in Eq. (J29j| with the 
expected values given in Eq. (|26|l are also shown in the table. 
The chi-square value per degree of freedom is xlof — 0.54. 



PCY-S1 O 



coupling const. 



value 



R 



(JVA) 

a s 

(NA) 

a v 
ANA) 
°s 

ANA) 
°V 

ai NA) 



-2.0305 x 10" 4 MeV" 2 0.79 



0.96 



1.6548 x 10" 4 MeV" 
2.2929 x 10" 9 MeV" 
-2.3872 x 10" 9 MeV 
-1.0603 x 10" 7 MeV 



1.37 



TABLE II. The parameter set PCY-S2 obtained by the omit- 
ting the tensor coupling. Xdof is 0.85. 



coupling const. 



value 



R 



ANA) 
a S 

ANA) 
°S 

ANA) 
°V 

ANA) 



-4.1595 x 10~ 5 MeV" 2 0.16 
1.3402 x 10" 5 MeV" 2 0.08 
1.3167 x 10" 9 MeV" 4 

-1.4018 x 10" 9 MeV" 4 




eA.j=;-i/2 — e Aj=i+i/2: when the A particle is put in 
the lowest s-orbital. Those obtained with PCY-S1 and 
PCY-S2 are shown in the lower panels of Figs. [T] and [2] 
respectively. For both the parameter sets, although the 
absolute values of A_E so are smaller by roughly a factor 
of 10 than those for nucleon, AE so alters its sign de- 
pending on the mass number. This does not happen in 
the meson exchange models (See Refs.[30l. l3ll|). and one 
may consider this inversion somewhat ill-favored. We 
mention, however, that at present there have been no ex- 
perimental data which exclude the possible inversion of 
the spin-orbit splitting in the medium and heavier mass 
region. 

If we exclude the spin-orbit splitting of the lp state of 
A in from the fitting, that is, if we fit only the energy 
ccntroid of each spin-orbit partner, we obtain Fig. [3] for 



TABLE III. The parameter set PCY-S3 obtained without fit- 
ting to the spin-orbit splitting in \ 6 0. Xdof i s 0.57. 



coupling const. 


value 


R 


(JVA) 
a S 


-2.0197 x 10" 4 MeV" 2 


0.79 


(NA) 

a v 


1.6449 x 10" 4 MeV" 2 


0.95 


ANA) 
°S 


2.3514 x 10" 9 MeV" 4 




ANA) 
°V 


-2.4993 x 10" 9 MeV" 4 




4 m) 


-4.0820 x 10" 9 MeV" 3 


0.05 
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FIG. 1. (Color online) Comparison between the experimental 
data and the calculated A binding energies Ba (upper panel) 
and spin-orbit splittings of A single-particle energies AE so 
(lower panel) obtained with the parameter set PCY-S1. The 
experimental data are taken from Refs. [4ll. [42l. |44| . 




0.12 0.14 0.16 0.16 



FIG. 2. (Color online) Same as Fig. [TJ but with the parameter 
set PCY-S2. 



the lambda binding energies and the spin-orbit splittings. 
The parameters for this set, PCY-S3, are summarized in 
Table UITl For this parameter set, the vector and scalar 
couplings of A to nucleon remain natural, but the ten- 
sor coupling is far smaller than the expected value in 
Eq. (|2"rj|) . Since there is no constraint on the value of 
spin-orbit splitting, this parameter set yields unaccept- 
ably large spin-orbit splitting, some of them stretching 
even beyond the experimental uncertainties (i.e., the up- 
per bounds). Again, a few of them have negative values. 

Lastly, we examine the role played by the derivative 
terms in the Lagrangian. In Ref. (45[, it was pointed 
out that only one derivative term is well constrained by 
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FIG. 3. (Color online) Same as Fig. [TJ but with the parameter 
set PCY-S3. 

TABLE IV. The parameter set PCY-S4 obtained by setting 
the vector derivative coupling, , to be zero, xlof i s 0-92. 



coupling const. value R 

a ( s NA) -1.8594 x 10" 4 MeV -2 0.72 
a { y A) 1.4981 x 1(T 4 MeW 2 0.87 

4 m) -1.9958 x 10" 10 MeV" 4 
6 ( V NA) 
a^ A) -5.5322 x 10" 8 MeV" 3 0.71 



the bulk nuclear observables, i.e., inclusion of a single 
derivative term is sufficient to obtain a good fit. Finelli 
et al. have shown that their model with only a scalar 
derivative coupling indeed reproduces well the data for 
normal nuclei [46j j and hypernuclei [34j . Following Finelli 
et al. [34], 5||, we construct another parameter set PCY- 
S4 by omitting the vector derivative term. The results 
are shown in Table IIVI and Fig. 2] One observes that 
the quality of the fit is as good as the other parameter 
sets. The agreement with the quark model prediction is 
also well, and the inversion of the spin-orbit partner is 
not seen for this force. Therefore, PCY-S4 provides an 
alternative parameter set to PCY-S1, where the main dif- 
ference between the two interactions is whether the spin- 
orbit splitting is normal (PCY-S4) or inverted (PCY-S1). 



IV. SUMMARY 

We have proposed a new relativistic point coupling 
model to describe single-A hypernuclei in the mean field 
approximation. This is a straightforward extension of the 
relativistic point coupling model for nucleons, which has 
a similar structure as the Skyrme interaction. To this 
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FIG. 4. (Color online) Same as Fig. [TJ but with the parameter 
set PCY-S4. 

end, we added effective contact N-A interactions, corre- 
sponding to the A- a and A-us couplings, to the model 
Lagrangian. In addition, we introduced the zero-range 
N-A tensor coupling as well to mimic the tensor cou- 
pling between A and uj meson, following the quark model 
suggestion. 

We fitted the coupling constants in the strange sec- 
tor to the experimental data of lambda binding energies. 
The four parameter sets, PCY-S1, PCY-S2, PCY-S3, and 
PCY-S4 were proposed, which well reproduce the exper- 
imental data through the whole mass region. The result- 
ing spin-orbit splittings in PCY-S1, PCY-S2, and PCY- 
S4 are smaller than that of nucleon by roughly a factor 
of 10 in their absolute values, although PCY-S3 yields 
too large spin-orbit splittings. For PCY-S1 and PCY-S2 
their signs are opposite to that of nucleon in some nuclei 
in the heavier region. On the other hand, for PCY-S4 
obtained without taking into account the vector deriva- 
tive term, the sign of the spin-orbit splitting is the same 
as that for nucleons. High-precision 7-ray experiments 
for A single-particle energies are awaited in order to see 
whether the spin-orbit splitting of heavy hypernuclei is 
normal or inverted. 

We have confirmed that the tensor coupling, which is 
ignored in the N-N interaction, is quite important to 
reproduce the small spin-orbit splittings of A particle. 
Without the tensor coupling, the scalar and the vector 
couplings of A to nucleon are forced to be unnaturally 
weak (PCY-S2). The tensor coupling suppresses the spin- 
orbit splittings, keeping the scalar and the vector cou- 
plings consistent with the naive quark counting. Those 
good consistency with the quark model found in our in- 
teraction can be a useful guide in further extending the 
point coupling model to multi-A or 3 hypernuclei. 

We conclude that the point coupling model is capable 
of describing single-A hypernuclei as well as normal nu- 
clei. The model can be an appropriate tool for relativistic 
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calculations for hypernuclei on three-dimensional mesh 
due to its numerical simplicity. There has been no rela- 
tivistic calculation performed on 3D mesh because of the 
variational difficulty Thus the primary future work is to 
develop an efficient calculation technique for relativistic 
calculations on 3D mesh that overcomes the "variational 
collapse" . A work in this direction is in progress. Further 
extensions of the point coupling model to multi-A and 3 
hypernuclei, and an introduction of explicit density de- 
pendences into the coupling constants are also interesting 
future works. 
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Appendix A: Center of Mass Energy 

In this Appendix, we give an explicit expression for 
the numerator of Eq. (|23|) . Although it has already been 
given in Ref. [471 ]. we have found a few typos in their 
expression for the relativistic case with spherical symme- 
try. Here we shall correct the typos and give the correct 
formula. 

With the spherical symmetry, single-particle (s.p.) 
wave functions are given as 



Vw (r) = 



(Al) 



a m 



where a is a shorthanded notation for {n ai , j Q }, 
while m — j z . Here, n a is the principal quantum num- 
ber, £^ is the orbital angular momentum of the upper 
component of the s.p. spinor, and j a is the total angu- 
lar momentum. The spherical spinor &tj m is defined by 
%m = T, m ' m "( im ' '\m"\jm)Y hn x^ m n , where X i m " is 
the spin wave function. The orbital angular momentum 
of the lower component is given by £^ = 2j a — ta ■ 
Following Ref. [47J , we compute the center of mass cor- 
rection in the non-relativistic approximation. The expec- 
tation value of the squared center of mass momentum 

)2 

CM 



(A2) 



fniu with a mean field many-body state reads 



a,0 



with 



77=±" 



and 



^<Vw|A|Vw) = (2j Q + 1) J2 i dr ^ 

KVw|vhW)i 2 

rn,m' 

= (2^ + 1X2^ + !)£(-) 



(v) 



d 2 [ 2 d g(g+l) 
dr 2 r dr r 2 



+2 



(v) 



(A3) 



77,77' 




(A4) 



T a 



r 



where w a is the occupation probability of the level a. 
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